< 

00 



Dh: 






o 



Solar Physics 
DOI: 10.1007/* 



Modification of Proton Velocity Distributions by 
Alfvenic Turbulence in the Solar Wind 



V. Pierrard^'^ • Y. Voitenko^ • 

IJ^ ' Received ; accepted 

' _ ' ^ © Springer •••• 

(N 

5_^ ■ Abstract In the present paper, the proton velocity distribution function (VDF) 

fS il in the solar wind is determined by solving numerically the kinetic evolution 

equation. We compare the results obtained when considering the effects of ex- 
ternal forces and Coulomb collisions with those obtained by adding effects of 
Alfven wave turbulence. We use Fokker-Planck diffusion terms due to Alfvenic 
turbulence, which take into account observed turbulence spectra and kinetic 
"^ . effects of finite proton gyroradius. Assuming a displaced Maxwellian for the 

proton VDF at the simulation boundary at 14 solar radii, we show that the 



QJ . turbulence leads to a fast (within several solar radii) development of the anti- 



sunward tail in the proton VDF. Our results provide a natural explanation for 

O . the nonthermal tails in the proton VDFs, which are often observed in-situ in the 

c/2 ' solar wind beyond 0.3 AU. 

c/j . 

Q , Keywords: Solar wind; Waves; Turbulence; Space Plasmas 

>>: 

f~i ■ 1. Introduction 

The solar wind is a low density plasma in which kinetic processes prevail. Kinetic 
models based on the solution of the Fokker-Planck equation have been developed 
^ . to study the steady state electron velocity distribution function (VDF) in the 

'^ I corona and at larger radial distances in the solar wi nd (SW) (|Lie-Svendsen et ai, 1997H 

^^ '. [Pierrard et al., 1999|[Lie-Svendsen and Leer, 2000[|Pierrard et ai, 2001[|Vocks and Mann, 2003H 

^^1 I 

(N 



Vocks and Mann, 2009). The test electrons were submitted to the influence of the 



external forces and to Coulomb collisions with background particles ( Pierrard et ai, 1999H 

|Lie-Svendsen and Leer, 20001 [Pierrard et al., 2001[ ). Such models were solved 

I^ . numerically and emphasized the effects of Coulomb collisions compared to the 

results of purely exospheric models that considered only the external forces 
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( jMaksimovic et ai, 1997| |Lainy et ai, 2003] ). They showed that Coulomb col- 
hsions have important effects on angular scattering (i.e. on the pitch angle 
distribution of the electrons) but do not modify their average density and mean 
temperatures radial distributions. Models including Coulomb collisions gave 
more realistic VDFs and a reduction of the temperature anisotropies compared 



to the purely cxosphcric approximation (Lemaire and Pierrard, 2001), because 
the Coulomb collisions isotropizc the electron VDF. Terms due to whistler wave 
turbulence have also been included in the electron kinetic equation and its in- 
fluence on the formation of nonthermal tails and on the velocity-space diffusion 



of electrons has been investigated (Pierrard et ai, 2011 1 



In the present paper, we study the evolution of the proton VDF in the solar 
wind. First, wc use a kinetic approach similar to the mentioned above papers, 
which accounts for the effects on the protons of the external forces and Coulomb 
collisions. After that wc include also the wave effects on the proton VDF. Recent 
in-situ measurements have revealed that the wave activity at the proton ki- 
netic scales is dominated by Alfven waves rather than whistlers ( |He et ai, 201 1[ 
|He et ai, 2012[podesta and Gary, 2011[ ). Among Alfven waves, the most power, 
about 80%, is in highly oblique kinetic Alfven waves (KAWs), and the rest 20% 
is in quasi-parallel ion-cyclotron Alfven mode ( |He et ai, 201 1[ ). Both the MHD- 
scalc and kinetic-scale Alfven waves posses power law spectral distributions and 
exhibit turbulence properties. We therefore use in present simulations a 1D-2D 
(ID spatial - 2D velocity-space) Fokker-Planck diffusion term due to Alfvenic 
turbulence derived by Voitenko and Pierrard (2012). 

Voitenko and Goossens (2006) suggested that the proton trapping in the 
parallel KAW potential can be responsible for the proton beams observed in the 
solar wind. Later on, at the Solar Wind 12 in June 2009, Pierrard and Voitenko 
have presented analytical estimations and numerical example, published in Pier- 
rard and Voitenko (2010). Numerical simulations by Osmane et at (2010) and 
Li et al. (2010) have demonstrated the proton beam formation by a coherent 
monochromatic KAW. Our present approach is based not on monochromatic 
KAWs but on a wide spectrum of KAWs with overlapping potentials, which 
results in the velocity-space proton diffusion rather than the proton trapping as 
in the mentioned above papers. 



2. Description of the Model 

The kinetic transport equation for the evolution of the velocity distribution 
function /(r, v, t) of the protons in the SW is given by: 

^^(^' ^' *) + (v.Vr)/(r, V, t) + (a. V.)/(r, v, t) = ('A + (A (1) 



dt \dt J ^ \dt ^ 

where r and v are respectively the position and velocity vectors of the particles, 
a is the acceleration due to external forces and t is the time. We are interested 
by the steady state solution of this equation. The right-hand-side represents the 
interactions between the particles. In the present paper, we consider the Coulomb 
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collisions with the Fokker-Planck term ( -^ J and wave-particle interactions with 

kinetic Alfven waves ( -j^ ) • 

In the case of the SW, the forces are the electric force ZeE, the gravitational 
force TOg and the Lorentz force resulting from the magnetic field distribution. 
The acceleration term for the protons is thus given by: 

/ZeE \ Ze, ^, , , 

a= +g + — vxB =a^r+aL 2 

\ m J m 

where Ze is the electric charge of the protons {Z = 1), m their mass, E the 
electric field, g the gravitational acceleration, B the magnetic field (assumed to 
decrease as r^^), a,.(r) is the radial component of the acceleration due to the 
electric and gravitational forces that are vertical and aL is the non radial term 
due to the Lorentz force. 

The simplest solar wind models consider only the effects of these external 
forces. In such models called exospheric, the interaction terms in Equation ([T]) 
are neglected so that an analytic solution of the equation can be obtained. 

Figure 1 shows the profiles of the VDF moments obtained with the exospheric 
solar wind model developed by Pierrard and Lemaire (1996) and improved by 
Maksimovic et al. (1997) and Lamy et al. (2003). In such models, truncated 
Maxwellian and Kappa VDF are assumed respectively for the protons and the 
electrons with given density and temperatures at the exobase level (here chosen 
to be at 1.5 R^ (solar radius)). The VDF is obtained at larger radial distance by 
solving the evolution equation, i.e. Equation ([1} without interaction right-hand- 
side terms. 

It can be seen (Figure Id) that the particles are accelerated at low radial 
distances due to the effects of the external forces. It is mainly from 1 to 20 
i?s that the bulk velocity of the particle increases and then it becomes almost 
constant at higher radial distances. This result is similar to that obtained in 
MHD solar wind models where the interactions between the particles are taken 
into account ( [Parker, 1958[ ). 

The exospheric model emphasizes that it is the electrostatic potential (see 
Figure lb) that increases the velocity of the particles. This potential ensures 
that the proton flux is everywhere equal to that of the electrons so that no net 
current is transported by the solar wind. 

In this model, the acceleration is reinforced by the presence of suprathermal 
electrons simulated by a Kappa VDF. They increase the electrostatic potential 
and thus accelerate the wind. On the contrary, the presence of suprathermal pro- 
tons has almost no influence on the electrostatic potential ( [Maksimovic et al., 1997[ ) .| 
That is why a Maxwellian is assumed at the reference level for the proton VDF 
at 1.5 i?s- The averaged values of the different moments (profiles of number 
density, bulk velocity, temperatures) are in good agreement with the solar wind 
observations, except the temperature anisotropics that are too high in the model 
and would be expected to produce plasma instabilities in the system. Such high 
anisotropics arc mainly due to the assumption of magnetic moment conservation 
for each particle species that is not realistic, even in low density plasmas. 
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Figure 1. Profiles of (a) number density, (b) electric potential, (c) flux, (d) bulk speed, (e) 
proton temperature Tp (upper solid line: Tm, bottom solid line: Tj_, dotted line: average T), (f) 
electron temperature (idem), (g) proton potential, and (h) heat flux (upper and lower curves 
for electrons and protons) from ro= 1.5 Rs up to 200 Rb obtained with a kinetic exospheric 
model of the solar wind with To = 1.5 X 10® K and Tp = 1.5 X 10® K for a Maxwellian proton 
VDF and a Kappa electron VDF with k=3. 



3. Numerical Model 



More sophisticated solar wind models including interaction terms (the right- 
hand-side of Equation (1)) need to be solved numerically. We develop in the 
present paper such a numerical model for the protons, based on similar assump- 
tions as those used for the model of SW electrons ( [Pierrard et ai, 1999[ ). 

We simplify the equation by using only three coordinates: the radial distance 
r, the velocity v and rj = cos 9 = Bq • k where 9 is the pitch angle between 
the velocity vector and the magnetic field direction assumed to be radial. The 
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left-hand-side of Equation (1) becomes: 

dt dr \ dv V di] J r drj 

As a first study, we will consider that the term {df/dt)c on the right-hand-side 
of Equation ([T|) represents the effects of the Coulomb collisions ( [Pierrard et ai, 2001[ ).| 
Then, in a next section of this paper, we adopt in addition the kinetic Alfven 
wave term [df /dt)p^. 



4. Spectral Method of Solution 

To obtain steady state solutions for the proton VDF, we use a specialized spectral 



method similar as that described in ( Pierrard, 2011 ) for the resolution of the 
Fokker-Planck equation (FPE) . 

The dimensionlcss velocity is defined by 

(4) 



2fcBTp(r) vth{r)' 

where wth is the proton thermal speed and k^ is the Boltzmann constant. 
The solution is expanded in terms of orthogonal polynomials: 

CJ-lN-lM-l \ 

^ ^ ^ aismPi{vi)Ss{x)L,n{z) j . (5) 

/=0 s=0 m=0 / 

where J, N , and M are integers whose value is adjusted to obtain the required 
numerical precision for the solution. We use Legendre polynomials Pi(rii) with 
respect to 77 = cos6', speed polynomials Ss{x) with respect to the normalized 
velocity x and displaced Legendre polynomials on the interval [0, c] with respect 
to the dimensionlcss altitude z. The aism coefficients arc the coefficients of the 
expansion to be determined. 

The boundary condition determines the value of some coefficients aism at the 
reference level. The equation is solved numerically to find the VDF f{z,x,r]) at 
the other radial distances. We choose 10 polynomials for each variable so that 
the results have a precision better than 10 % with reasonably short CPU times. 
We have checked that the results do not change significantly when J, N 01 M 
is increased. 

The advantage of this method is that in a discrete ordinate basis, the deriva- 
tives of any continuous function f{y) can be approximated by the following 
expansion: 

flf\ ^^^ 

^ ' V=Vi j=Q 

where Dij are the matrix elements of the derivative operator in the polynomial 
basis. 
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The integrals to calculate the moments of the protons are also easily per- 
formed by numerical quadrature taking into account the polynomial expansion 
( [Pierrard and Lemaire, 1996[ ). The moments are calculated from the VDF found 
by solving the FPE. The density is given by: 



n{r) = / /(r,u,r;) dv. (7) 

The bulk velocity is calculated by: 

u(0 = ( / fir,v,r])vdv ] /n. (8) 



The temperature corresponds to a second order moment of the VDF: 

Tir)^^mir)+2T^ir)) (9) 

where 

T||(r) = -^ ^ " , ; (10) 

„ / N m f v^, f(r, V, n) dv , ^ 

^-^W = ^ ^i^ ) ' • 11) 

2fcB7i(r) 



5. The Effects of Coulomb Collisions 

The effects of the binary Coulomb collisions are taken into account by using as 
the right-hand-side term in Equation ([1]) the Fokker-Planck collision operator 
appropriate when large-angle deflections can be neglected (Hinton, 1983): 



df\ d 



dt J ^ 9v 



A/(r,v,t)-i|-.(D/(r,v,t)) 



(12) 



where A is the dynamic friction vector 
A = -47r 



ZlZJ e^lnA , rrir 

171/3 



X l^-'M-')^^ (13) 



and D is the velocity diffusion tensor 



D = Att^^^ (14) 



x/^v'Mv')(^ 



(v-vO(v-V) 
{v — w')-^ 
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In Equations ([T^ and p^ . In A is the usual Coulomb logarithm containing 
the Debye screening effect: InA « 24. The index a corresponds to the test 
protons and /3 to the background particles, m the mass and Ze the charge of the 
particles, ffi{v') is the velocity distribution function of the background particles 
assumed to be a displaced Maxwellian distribution since this VDF is solution of 
the Fokker-Planck operator. 

We assume boundary conditions at a reference level here chosen to be 14 R^. 
It could be chosen at lower radial distances as well, but this distance corresponds 
to the region where kinetic Alfven turbulence is supposed to have effects, so this 
altitude is especially appropriate for comparisons when this term included (see 
the next section). We assume an isotropic displaced Maxwellian for the protons 
VDF at this reference distance since this VDF is solution of the Fokker-Planck 
Coulomb collision term. The bulk velocity at this radial distance is chosen to be 
470 km s^^, as obtained with the exospheric model. 

Figure 2 illustrates the VDFs (normalized to the thermal velocity and centered 
on the averaged bulk velocity) obtained at different radial distances by solving 
the FPE including the Coulomb collision term. The VDF is always presented in 
the reference frame centered at its bulk velocity. The dotted circle corresponds 
to w/wth = \/372- 

The VDF remains close to a displaced Maxwellian at higher radial distances 
with a bulk velocity slightly increasing as illustrated on Figure 3b. Figure 3a 
illustrates the profile of the decreasing normalized number density and Figure 
3c the decreasing temperature. Compared to the exospheric model, the effects 
of Coulomb collisions (obtained with the inclusion of the Fokker-Planck colli- 
sion term and with the isotropic boundary condition) are mainly to reduce the 
temperature anisotropy of the particles. 

The results depend on the boundary condition at the reference level. Since we 
start with an isotropic displaced Maxwellian at the reference level chosen to be 14 
Rs, the VDF remains close to a displaced Maxwellian at higher radial distances 
even if the Coulomb collisions become negligible at large radial distances. The 
bulk velocity increases with the radial distance due to the acceleration term and 
the Lorentz force slightly alignes the VDF to the magnetic field direction with 
the increasing radial distances. 



6. Fokker-Planck DifTusion Terms due to Alfvenic Turbulence 

In this section we include diffusion terms due to Alfvenic turbulence. The Alfve- 
nic turbulence evolves towards small scales across the background magnetic field 
and at proton kinetic scales transforms into the turbulence of KAWs. The pre- 
liminary principles of KAWs on the formation of the proton beams have been 
presented by IPierrard and Voitenkol (|2010p . The focus in this previous paper 
was on the proton trapping and acceleration in the electric potential carried 
by intermittent isolated KAWs. Here we use a model Alfvenic spectrum that 
corresponds to the magnetic power spectra measured in the solar wind. 

The Alfvenic turbulence and related terms describing the velocity-space pro- 
ton diffusion in the solar wind conditions have been presented in detail by 
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a) 
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a) 
o 
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Figure 2. VDFs obtained at diflFcrcnt radial distances assuming a displaced Maxwellian as 
boundary condition at 14 Rb by solving the FPE including the Coulomb collision term. 



Voitenko and Pierrard (2012). We use here the final expressions from this paper, 
which are in terms of quantities measured in-situ in the solar wind. The Fokker- 
Planck term due to Alfvenic turbulence, which is essentially the KAW turbulence 
at small scales, is given by 



dfp 
dt 



dv V dij J \ dv 



I -if d 



2—D —f+ —l—^D —f 
dv dv dv V Ot] 

l-if d ^ d ^ 1 - if d 1 - 

+ TT'^^^Trf + a 

V df] dv V drj V 



d-q 



drj 
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Figure 3. Radial profiles of number density, bulk velocity, and proton temperature profiles 
obtained by solving the FPE including the Coulomb collisions assuming a displaced Maxwellian 
as boundary condition at 14 Rs. 



The KAW diffusion coefficient 



Da = aTT^flpYsVA 



[ S ) \Bi\- 



{{vv/VAf + l) (vv/Va) 



B2 



for 1 < -qv/VA < ■\/l + 2/1^^ and Da — otherwise, where /z — k±pp is the 
perpendicular wavenumber k± normahzed by the proton gyroradius pp. Other 
notations are as follows: 51p is the proton cyclotron frequency. Va is the Alfven 
velocity: 



VA(cm s-^) = 2.18 X 10 



11 B (gauss) 



•y/n(cni •^) 



(15) 
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Vs = -v/fce Tell /nip is the ion-acoustic velocity, and the anisotropy factor a ~ 0.3. 
All the plasma parameters (background number density no, parallel electron 
temperature T^w, and background magnetic field Bq) are functions of r. 

The above expression for Da is obtained from Equation (7) by Voitenko and 
Pierrard (2012) using a power-law turbulence spectrum |i?^sc| oc |i?i| /i"''^ 
with the following model for the variable index: S± = 1.6 + 1.1/i^ (l 4- /i^) 
This S± has the low- and high-/i asymptotes 1.6 and 2.7, as is measured in the 
solar wind (Alexandrova et at, 2009). We also averaged Da in the cross-field 
plane and put T^/Tp = 1 in the KAW dispersion. 

A low-beta two-fluid dispersion equation is adopted for KAWs, 

" ^ (16) 



k.VA 




which is a good approximation for the kinetic dispersion by Hasegawa and 
Chen (1976) where a possible temperature anisotropy is taken into account. 
A comparison among several KAW models is given by Voitenko and Pierrard 
(2012). 



7. Formulation of the Boundary- value Problem at the Solar Wind 
Base 

We are interested in the quasi-stationary boundary-value problem, where the 
MHD Alfvcn waves (AWs) are injected along Bq at r = tq = 14 i?g. We assume 
-6(200 i?s) = 4 nT = 4 x 10~^ gauss and decreasing as r^^ to determine B 
at 14 Rs- Then d/dt means v^d/dr in the Fokker-Planck Equation ([1]). The 
nonlinear interaction of MHD AWs create MHD turbulent cascade towards 
smaller scales and turbulent spectra. It is found that the turbulent cascades are 
anisotropic and proceed mainly towards small cross-field wavelengths, where the 
waves are essentially kinetic Alfven waves undergoing a wave-particle Cherenkov 
interaction. 

As the plasma parameters, magnetic field -Bq, and Alfven velocity Va vary 
with distance, we expect that the upward propagating KAW turbulence sweeps 
the anti-sunward part of the particle distribution functions making them flatter 
or even producing a separate beam. Such distortions of the velocity distributions 



arc often measured at the distances > 0.3 AU ( Marsch et al, 1982 ), which can 
be signatures of the velocity diffusion due to the KAW turbulence. 

We assume that, contrary to the distant solar wind, the Alfvcnic turbulence 
in the vicinity of the solar wind base does not fill all the wavenumber range but 
evolves rapidly. We will model this by introducing a distance-dependent high- 
wavenumber cutoff /ic for the spectrum of such "young", not fully developed 
turbulence. The turbulent spectrum in the "young" turbulence does not spread 
beyond /i = /ic, but fic itself increases with heliocentric distance. Because of the 
absence of the observational data on the cutoff /ic within 0.3 AU, we model it as 

3 
yUc = Mc (r) = /fcO — 
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where //cO is the cutoff value at the solar wind base r ~ tq. We assume fico = 0.1 
at the boundary, such that the turbulent spectrum does not reach the proton 
kinetic length scales w prmp- The proton gyro-scale /ic ~ 1 is reached by the 
turbulent cascade at the distances sa 2ro. 

In the near-Earth solar wind at i? = 1 AU the normalized spectral density 
1^1 1 /Bq « 10"'^ (Alexandrova et ai, 2010). Since the radial dependence of 
1^1 1 / Bq is not known, we apply a self-similar projection from 1 AU to tq = 14 
i?s by keeping constant the ratio of the total turbulence energy to the spectral 
density at yu = 1. Then we use a model profile 



l^^l^ 0.05 



B'o 



0.68 



(17) 



describing a relatively slow decrease of \Bx\ /Bq with r starting from l-Bil"" /Bq 
w 1.6xl0~2atr = ro = 14i?s,to \Bif /B'^ w4xl0-3at 1 AU. Alexandrova et 
al. (2010) reported values > 10"'^ at 1 AU. This profile fits the turbulence energy 
in the time domains from 3 to 25 minutes, obtained by Bavassano et al. (1982) 
from Helios measurements at three radial distances between 0.3 and 1 AU. Other 
attempts to model the development of solar wind turbulence have adopted other 



turbulence profiles ( Cranmer and van Ballegooijen, 2003 Cranmer et al, 2007 1 , 



which imply qualitatively similar turbulence behavior above 15 i?s- 

The spectral power density of magnetic fluctuations at any scale fj, is W^ = 
B^/^, where B^ is the magnetic amplitude at that scale. The total energy of the 
turbulence is given by the integral over all ^ from the lowest cutoff at electron 
scales to the largest (driving) scale at the break in the spectral index between 
-1 and -5/3. Most of this energy is carried by the large-scale MHD waves that 
compose the MHD inertial range. Formally the energy in the steep kinetic spectra 
of KAWs immediately available for particles is much less, of the order of q~^^^ 
(here q = k±s/k±i is the driving-to-kinetic scale ratio, or "width" of the inertial 
range). The wave propagation angle at proton kinetic scales is close to 7r/2. Our 
Bi is an average fluctuating amplitude at a particular break scale ^ = ^i (we 
assumed fj.i=l for simplicity ). It is therefore much smaller than the total energy 
of Alfvcn waves at all scales. 



8. Modification of Proton VDF by the Alfvenic Turbulence 

It seems that the observed turbulence spectra do not evolve significantly beyond 
0.3 AU, which means they are formed much closer to the Sun. Since the kinetic- 
scale Alfvenic turbulence and its effects on protons are our primary interest, it 
is convenient to set the simulation boundary at the radial distance where the 
turbulence cascade approaches kinetic scales. As a reference level we chose here 
14 Rs where we allow the turbulence cascade to spread in the high-wavenumber 
range. Again, the solar wind is already formed in this region and the turbulence 
terms can be included in consideration without interfering much with the SW ac- 
celeration processes. We assume a displaced Maxwellian for the protons VDF at 
this reference altitude since this VDF is solution of the Fokker-Planck Coulomb 
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collision term. The bulk velocity at this radial distance is chosen to be 470 km 
s~^. This choice is conventional and the boundary could be at lower or larger 
radial distances, which however could not change simulation results except for 
adjusting them to a new distance. 

Figures 4-6 present results of the simulations including the turbulence terms 
for the radial distances starting from the simulation boundary at 14 Rg. At the 
boundary, the same displaced MaxwcUian for protons was used as in the previous 
simulations. The results presented in Figure 4 show the proton VDF evolution 
with the radial distance. As compared to the previous simulation without turbu- 
lence, where the proton VDF remains a displaced Maxwellian at all distances, the 
presence of the Alfvenic turbulence in the solar wind leads to a fast development 
of nonthermal tails in the proton VDFs. The tails are already noticeable at the 
distances about one solar radius from the boundary, and they increase fast with 
the radial distance and become pronounced after 2-3 solar radii. 

The radial behavior of the moments of the proton VDF are shown in the top 
panels of Figure 5. As is seen from Figure 5, the profiles of the bulk velocity 
and number density are also modified by the turbulence. The bulk velocity is 
increased in the regions where the tails are formed, and it is stabilized at higher 
altitudes. The number density profile decreases slower than without turbulence. 
The number density profile decreases less fast when the K AW amplitude Bi /Bq 
is chosen to be smaller. The bottom panels of Figure 5 show the profiles of the 
Alfven velocity and of the plasma beta function: 



Tr4 \ Vl 



^Hi+5trj^ ™ 



where Vxp = \/kBTp/mp. The most favorable conditions for the tail generation 
occur in the regions where the proton thermal and Alfven velocities are about 
the same, i.e. /3 « 1. 

Figure 6 illustrates the VDF obtained at 30 i?s assuming a displaced Maxwel- 
lian as boundary condition at 14 Rs by the evolution equation including the 
Alfvenic turbulence terms. The top panel shows the parallel (solid line) and per- 
pendicular (dotted line) cross section of the proton VDF. It can be seen that an 
extended nonthermal tail is created in the direction parallel to the interplanetary 
magnetic field. The Coulomb coUisional diffusion appeared to be efficient enough 
to smooth out a strong gradient at the turbulent diffusion boundary T4 w Va and 
to prevent formation of the quasilinear plateau aiVz ^ Va. The core distribution 
seems not to be affected. The bottom panel illustrates the isocontours of the 
proton VDF in the velocity plane centered on the bulk velocity. 

Particles VDF measured in-situ in the solar wind show many non-thermal 
features, like particle beams, thermal anisotropics, heat fluxes, etc., which are 
possible to explain with kinetic modeling. In particular, contour plots of the 
proton VDF measured around 0.3 AU by Helios showed the nonthermal tails 
and sometimes also separate bumps on tails (proton beams) formed at low radial 
distances ( |Marsch et ai, 1982| [Marsch, 2006] ). 



Our simulation results demonstrate that the kinetic-scale Alfvenic turbulence 
can create the nonthermal tails in the proton VDFs in the direction of the local 
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Figure 4. Isoeontours of the proton VDF obtained at diflFercnt radial distances assuming 
a displaced MaxwcUian as boundary condition at 14 Rb- The evolution equation was solved 
numerically and included the Fokker-Planck diffusion terms due to turbulence. 



magnetic field. The temperature anisotropy of the core is not obtained with this 
term, but could be by including the 2D diffusion term associated to KAW. 



9. Conclusions 

An essential difference between MHD and kinetic modeling is the ability of the 
latter to capture processes of the coUisionless wavc-particlc interactions. Kinetic 
effects arc of particular importance in weakly and mildly collisional media, like 
solar and stellar coronae and winds, planetary magnetospheres, and interstellar 
medium. In particular, permanently observed kinetic-scale Alfvenic turbulence 
need to be incorporated in kinetic solar wind modeling. 

We investigated numerically the influence of the kinetic-scale Alfvenic tur- 
bulence on the evolution of the proton distribution function and the principal 
moments in the SW at low radial distances. The effects of the Coulomb collisions 
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Figure 5. Radial profiles of (a) number density, (b) bulk velocity, (c) the Alfven velocity, and 
(d) the plasma beta function. Panels (a) and (b) are obtained by integrating the proton VDF 
weighted by 1 and Vz, respectively. Kinetic Alfven turbulence has been taken into account in 
the evolution equation assuming a displaced Maxwellian as boundary condition at 14 R^. 



and external forces are compared to those included also turbulence. Observed 
properties of Alfvenic turbulence are taken into account, including spectral 
breaks, variable spectral slopes, and radial dependence of the turbulence level. 
We found significantly different proton VDFs with and without turbulence. In 
particular, the obtained simulation results demonstrate that the Alfvenic turbu- 
lence with characteristics typical for the solar wind is very efhcient in producing 
the nonthermal tails in the proton VDFs in the direction of the background 
magnetic field. We observed in our simulations a fast generation of the nonther- 
mal proton tails at the heliocentric distances less than 0.3 AU. The tails do not 
appear without turbulence terms. This provides a natural explanation for the 
proton VDF with tails routinely observed in situ at 0.3-1 AU. 

The turbulence effects on the number density and bulk velocity are also 
investigated. It is shown that the turbulence has a non negligible influence on 
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Figure 6. VDF obtained at 30 -Rs assuming a displaced Maxwellian as the boundary condition 
at 14 Rb by solving the evolution equation including the Fokker-Planck diffusion term due to 
kinetic Alfvcn turbulence. Top panel: Parallel (solid line) and perpendicular (dotted line) cross 
section of the proton VDF. Bottom panel; Isocontours of the proton VDF in the velocity plane 
centered on the bulk velocity. 



the bulk velocity and increases it in the formation region of the proton tails. 
On the contrary, the particle number density decreases when the turbulence is 
accounted for. 

In this first kinetic simulation we considered a simplest case of the turbulence 
spectrum with a single spectral kink, which did not allow us to study the influ- 
ence of steeper spectra with indices Sj_ cc 3 — A often observed in the solar wind. 
Again, we used a simplest ID version of the diffusion coefficient reduced by the 
averaging in the cross-field velocity plane. 

The bump-in-the tail VDFs of protons (proton beams) , observed in the solar 
wind, are not easy to reproduce by the stationary Fokker-Planck diffusion. As 
is explained in our accompanying paper (Voitenko and Pierrard, 2012), a time- 
varying diffusion can lead to the beam formation by the time-of-flight effects. 
The beams can be also produced by the proton trapping and acceleration in 
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the parallel electric potential of the intermittent KAW pulses, as is explained by 
Pierrard and Voitenko (2010) and demonstrated numerically by Li et at (2010). 
Simulation of this process requires a more complex modeling of the solar wind 
turbulence including intermittence. This work is in progress. 
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